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The extensive study of nuclear multiframentation for the last two decades 
has been strongly stimulated by the idea that this process is related to a 
liquid-gas phase transition. One of the first nuclear models, suggested by 
N. Bohr, K. Weizsacker and Ya.I. Frenkel 65 years ago, is the liquid-drop 
model, which is still alive. A liquid-gas phase transition in nuclear matter 
was predicted much later [1-3] on the basis of the similarity between van 
der Waals and nucleon-nucleon interactions [2j. The equations of state for 
the two cases are similar. A point of particular interest is the existence of a 
spinodal region at reduced densities characterized by a phase instability. One 
can imagine that a hot nucleus expands due to thermal pressure and enters 
into this unstable region. Due to density fluctuations, a homogeneous system 
converts into a mixed phase, consisting of droplets (IMF's, 2 < Z < 20) 
surrounded by nuclear gas. In fact, the final state of this transition is a 
nuclear fog [3] ■ The neutrons fly away with energies corresponding to the 
system temperature (6 — 7 MeV), while the charged particles are additionally 
accelerated by the Coulomb field of the system. The disintegration time is 
determined by the time scale of the density fluctuations and is very short. 
Indeed, it was measured in number of papers that the IMF's emission time 
is less than 100 fm/ c. This is the scenario of nuclear multifragmentation as a 
spinodal decomposition, considered in a number of papers (see, for example, 
[4-12], and review papers [13-14]). The spinodal decomposition is, in fact, 
the liquid-fog phase transition in a nuclear system. 

An important parameter of this scenario is the critical temperature for 
the nuclear liquid-gas phase transition T c at which the isotherm in the phase 
diagram has an inflection point. The surface tension vanishes at T c , and only 
the gas phase is possible above this temperature. There are many calcula- 
tions of T c for finite nuclei. In Ref. [Tl 121 [TBI I16j . for example, it is done by 
using a Skyrme interaction and the thermal Hartree-Fock theory. The values 
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of T c were found to be 10 — 20 MeV, depending on the Skyrme interaction 
parameters and the details of the model. Experimental estimations of the 
critical temperature for the finite nuclei have been done in several papers. 

The main source of experimental information for T c is the fragment yield, 
but the procedures to extract T c are heavily debated. In some statistical 
models of nuclear multifragmentation the shape of the IMF charge distri- 
bution y(Z) is sensitive to the ratio T/T c . The charge distribution is well 
described by the power law Y(Z) ~ Z~ r for a wide range of colliding sys- 
tems |17j . In earlier studies on multifragmentation [3J I18j the power-law 
behavior of the IMF yield was interpreted as an indication of the proximity 
of the excited system to the critical point for the liquid-gas phase transi- 
tion. This was stimulated by the application of Fisher's classical droplet 
model which predicted a pure power-law droplet-size distribution with 
the minimal value of r = 2 — 3 at the critical point. 

In Ref. ^H] Hirsch et al. estimate T c to be ~ 5 MeV simply from the 
fact that the mass distribution is well described by a power law for IMF's 
produced in the collision of p (80 — 350 GeV) with Kr and Xe targets. In 
fact, the fragment mass distribution is not exactly described by the power 
law, therefore it was suggested the use of the term T app , an apparent expo- 
nent, to stress that the exact power-law description takes place only at the 
critical temperature. In paper [2U] the experimental data were gathered for 
different colliding systems to get the temperature dependence of T app . As a 
temperature, the inverse slope of the fragment energy spectra was taken in 
the range of the high-energy tail. The minimal value of T app was obtained at 
T = 11 — 12 MeV, which was claimed as T c . The later data smeared out this 
minimum. Moreover, it became clear that the "slope" temperature does not 
coincide with the thermodynamical one which is several times smaller. 

A sophisticated use of Fisher's droplet model for the estimation of T c has 
been recently made by Elliott, Moretto et al. |21|I22|. The model was modi- 



3 



fied by including the Coulomb energy release when a particle moves from the 
liquid to the vapor. The data from the Indiana Silicon Sphere Collaboration 
for 7r (8 GeV/c) + Au collisions were analyzed (2J. The extracted critical 
temperature is T c = 6.7 ±0.2 MeV. In a recent paper [22j the same analysis 
technique is applied to the data for the multifragmentation in collisions of 
Au, La, Kr (at 1.0 GeV per nucleon) with a carbon target. The extracted 
values of T c are 7.6 ± 0.2, 7.8 ± 0.2 and 8.1 ± 0.2 MeV respectively. 

There is only one paper in which T c is estimated by using data other than 
the fragmentation ones. In Ref. |23| it is done by the analysis of the temper- 
ature dependence of the fission probability for 4 He + 184 W collisions [21] ■ It 
was concluded that T c > 10 MeV in contrast to the result of Ref. |21[ I22j. 

It should be noted that in some papers the term "critical temperature" 
is used in another meaning than given above. In Ref. [2E] multifragmenta- 
tion in Au + Au collisions at 35 A MeV was analyzed with the so-called 
Campi plots to prove that the phase transition takes place in the spinodal 
region. The characteristic temperature for that process was denoted as 
T cr n and found to be equal to 6.0 ± 0.4 MeV. In the recent paper [2E] the 
bond percolation model is used to interpret 10.2 GeV/c p + Au multifrag- 
mentation data. The critical value of the percolation parameter p c = 0.65 
was found from the analysis of the IMF charge distribution. The corre- 
sponding "critical temperature" of 8.3 =L 0.2 MeV is estimated by using the 
model relation between the percolation control parameter u p" and the ex- 
citation energy. The more appropriate term for this particular temperature 
is "break-up" or "crack" temperature, as suggested in Ref. [22 • This tem- 
perature corresponds to onset of the fragmentation of the nucleus entering 
the phase coexistence region. The low-multiplicity channels dominate dur- 
ing the onset of multifragmentation characterized by a U-shaped fragment 
mass distribution. As shown by means of the statistical multifragmentation 
model (SMM) |27| I28j. the average hot fragment multiplicity is M=3-5 at an 
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excitation energy around 4 MeV/nucleon, and the probability of the com- 
pound nucleus channel is still considerable. Exactly at these relatively low 
excitation energies the experimenters observed the critical phenomena (see, 
for example [HI [23 E3 ) • 

Having in mind the shortcomings of Fisher's model |311 I32| . we have 
made an attempt to estimate the critical temperature in the framework of 
SMM. It describes well the properties of the thermal fragmentation of target 
spectators produced in collisions by light relativistic ions. As an example, 
Fig. 1 (top) shows the fragment charge distributions measured by the FASA 
Collaboration for collisions of p (8.1 GeV), 4 He (4 and 14.6 GeV) and 12 C 
(22.4 GeV) with Au targets ^2j along with the calculated charge distribu- 
tions. The mechanism for the reactions of light relativistic projectiles is 
usually divided into two stages. The first is a fast energy-depositing stage, 
during which very energetic light particles are emitted and a nuclear remnant 
is excited. We use the intranuclear cascade model (INC) [33] for describing 
the first stage. The second stage is described by SMM, which considers the 
multibody decay of a hot and expanded nucleus. But such a two-stage ap- 
proach fails to explain the observed IMF multiplicities. An expansion stage 
(Exp) is inserted between the two parts of the calculation. The excitation 
energies and the residual masses are then fine tuned on an event-by-event 
basis |12j to get agreement with the measured IMF multiplicities. The lines 
in Fig. 1 (top) give the charge distributions calculated in the framework of 
this combined model assuming T c = 18 MeV. The agreement between the 
data and the model prediction is very good. The bottom panel of Fig. 1 
shows the power-law fit of the distributions with the r parameter given in 
the insert as a function of the beam energy. The corresponding thermal 
excitation energy range is 3 — 6 MeV/nucleon. The power law parameter 
exhibits the so-called "critical behavior" showing a minimum at an exci- 
tation energy corresponding to a temperature three times lower than the 
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assumed T c . A conventional explanation of the occurrence of a minimum is 
given in Ref. QUE]. 

The charge yield depends on the contribution of the surface free energy 



This equation was obtained in Ref. [31] devoted to the theoretical study 
of thermodynamical properties of a plane interface between two phases of 
nuclear matter (liquid and gas) in equilibrium. The corresponding calcu- 
lations were performed with the Skyrme interaction. The phase diagram 
generated by the SMM model (using eq. (1)) is discussed in detail in |35j . 
This parametrization is successfully used by the SMM for describing the 
multifragment decay of hot finite nuclei. In particular the SMM describes 
the experimental critical behaviour of fragments and scaling in multifrag- 
mentation |25l 1291 I3U| with the standard T c = 18 MeV. This scaling was 
taken as a starting point of the analyses |21[ I22| also. 

The present calculations are performed for p (8.1 GeV) + Au collisions 
with T c as a free parameter. For all values of T c the calculations with the 
INC+Exp+SMM model have been properly adjusted ^2] to get the mean 
IMF multiplicity close to the measured one. Figure 2 (left) shows the com- 
parison of the measured fragment charge distribution with the model predic- 
tions for T c = 7, 11 and 18 MeV. The statistical errors of the measurements 
do not exceed the size of the dots. The calculations are close to the data 
for T c = 18 MeV. The estimated mean temperature of the fragmenting sys- 
tem is around 6 MeV, the mean charge and mass numbers are 67 and 158 
respectively. The theoretical curves deviate from the data with decreasing 
T c . The right panel gives the results of the power-law fits for the data and 
model calculations (in the range Z=3 — 11). 



of the fragments (a s (T)A 2 ^ 3 ) to the entropy of a given final state of the 
partition. The following expression is used in the SMM for a s (T): 




(1) 
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The final results are shown in Fig. 3. The measured power-law exponent 
is given as a band with a width determined by the statistical error. The size 
of the symbols for the calculated values of T app is of the order of the error 
bar. The model predicted values of the power-law exponent are significantly 
smaller than the measured one for the range of T c < 13 MeV. From the best 
fit of the data and calculations one concludes that T c = 20 ± 3 MeV (90% 



Figure 3 shows also the results of the calculations with a s (T) linearly 
dependent on T/T c [STIISS]: 



The calculated values of r app in this case are remarkably lower than 
the measured one for any value of the critical temperature used (up to 
T c = 24 MeV). 

To conclude, the IMF charge distribution for p + Au collisions at 8.1 GeV 
has been analyzed within the Statistical Multifragmentation Model with 
T c (at which surface tension vanishes) as a free parameter. The value 
T c = 20 ± 3 MeV (90% CL) obtained from the best fit to the data is con- 
sidered as an effective value of the critical temperature averaged over all 
the fragments produced in the collision. This value is significantly larger 
than those found in Ref. |21l I22j by the analysis of the multifragmentation 
data in terms of Fisher's droplet formalism. A surprisingly large range of T c 
values in different publications indicates on the severe model dependence of 
the results. Although our value for T c is model dependent, as is any other 
estimate of the critical temperature, the analysis presented here provides 
strong support for a value of T c > 15 MeV. Another conclusion which can 
be drawn from this work is that the properties of individual hot fragments, 
in particular, their surface energies, can be obtained from the experimental 
data, and they are extremely important for identification of the phase tran- 



CL). 
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sition. This puts additional constrains on models used for description of the 
phase transitions in nuclear systems. 
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Figure 1. Karnaukhov, PRC. 
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Figure 3. Karnaukhov, PRC. 
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Figure Captions 



Fig. 1. Fragment charge distributions for p+Au at 8.1 GeV, 4 He+Au at 4 
GeV, 4 He+Au at 14.6 GeV and 12 C+Au at 22.4 GeV. The lines (top) are 
calculated and normalized at Z=3. The power-law fits are shown on the bot- 
tom panel with r parameters given in the insert as a function of beam energy. 

Fig. 2. Fragment charge distribution for p+Au at 8.1 GeV (dots). The 
lines (left side) are calculated assuming T c = 18 MeV (solid), 11 MeV (dot- 
ted) and 7 MeV (dashed lines). The power law fits are presented on the 
right panel. 

Fig. 3. The power-law exponent for p(8.1 GeV)+Au. The band corre- 
sponds to the measured value and its error bar. The symbols are obtained 
by the power-law fits of IMF charge distributions calculated assuming dif- 
ferent values of T c and different parameterizations of the surface tension: 
squares are for eq. (1), solid circles are for eq. (2). 
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